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£H | Abstract 

We consider Lagrangian systems in the limit of infinitely many par- 
ticles. It is shown that the corresponding discrete action functionals 
Gamma-converge to a continuum action functional acting on probabil- 
ity measures of particle trajectories. Also the convergence of stationary 
points of the action is established. Minimizers of the limiting functional 
and, more generally, limiting distributions of stationary points are investi- 
| gated and shown to be concentrated on orbits of the Euler-Lagrange flow. 

^ (-h We also consider time discretized systems. These results in particular 

provide a convergence analysis for optimal transportation meshfree meth- 
ods for the approximation of particle flows by finite discrete Lagrangian 
dynamics. 
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> ! 1 Introduction 

in 

In classical Lagrangian mechanics a system of particles is described by an action 
| functional on the particle trajectories. According to Hamilton's principle the 

dynamics of the system is given by stationary points of this functional. In 
continuum mechanics, on the other hand, a medium such fluid is 

described by its mass density and fields for physical observables whose dynamics 
are governed by partial differential equations such as the velocity obeying, e.g., 
Euler's equation of motion. While at microscopic molecular length scales a 
fluid cannot be viewed as a homogeneous medium, these fields are assumed to 
describe material properties in a suitably mesoscopically averaged sense. The 
natural question therefore arises, if these different models can be related to each 
other. 

In order to address this question, one is naturally led to extend the set-up 
of Lagrangian mechanics to infinite dimensional systems and to devise action 
functionals acting on mass densities or, more generally and after normalization, 
probability measures that describe the mass distribution of the system. If with 
the help of a single particle Lagrangian a 'Lagrangian cost function' is defined 
as the minimal value of the action necessary to move a particle from one point in 
space to another, a variational framework for such systems can be given within 
the theory of optimal transportation. (See, e.g., [12] for a recent account on 
optimal transportation theory with Lagrangian costs.) 

A connection between mass transportation problems and fluid mechanics 
has been derived in the seminal work [21 [3] of Benamou and Brenier (see also 
[llj). who show that the Wasserstein distance between two probability densi- 
ties can be expressed as a minimal action value with Lagrangian cost, when 
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minimized over velocity fields constrained to satisfy the continuity equation 
ensuring conservation of mass. Departing from this relation, Li, Habbal and 
Ortiz have constructed an efficient Lagrangian meshfree approximation scheme 
for fluid flows by material point sampling and solving for discrete minimizers of 
an approximating finite dimensional Lagrangian system, see [8]. An interesting 
aspect of this approach is that, due to their variational structure, the discrete 
approximating systems can be formulated within the theory of variational in- 
tegrators and consequently possess good conservation properties of physically 
conserved quantities. See, e.g., [9] and the references therein for a general in- 
troduction to the theory of variational integrators and [10] for a convergence 
analysis on manifolds. As discussed in [8], conforming fields may be interpo- 
lated efficiently with max-ent shape functions as developed in [I], for which 
convergence has been proved in the recent article [6]. 

The main aim of the present work is to provide a rigorous variational analysis 
of the infinite particle limit in Lagrangian mechanics. Such an analysis appears 
to be interesting from at least three different perspectives: 

Firstly our results supply a rigorous convergence analysis of the optimal 
transportation meshfree method constructed in [8]. We refer to [8] for an ex- 
tensive description of this method, its comparison with alternative numerical 
integration schemes and numerical experiments. 

Secondly our approximation scheme also provides a novel method to infer 
characteristic properties of optimal transport maps and their displacement in- 
terpolation with Lagrangian costs. In particular we show that with the help of 
our discrete-to-continuum approximation scheme through T-convergence it is 
possible to re-derive a result of Bernard and Buffoni [1] characterizing dynam- 
ical optimal transference plans. 

Thirdly the problem is of some independent interest as it provides an effec- 
tive theory derived by a rigorous variational discrete-to-continuum r-limit for 
a time dependent problem. Such problems have received a lot of attention over 
the last years in particular for static problems, see, e.g., the survey article [5] 
by Blanc, LeBris and Lions. It has to be noted, however, that the much more 
interesting problem for interacting particles is beyond the scope of this paper 
and deferred to future work. 

Having reviewed some basic material on Lagrangian systems for single and 
many particle systems in Section [2j in Section 13.11 we state and prove the main 
variational convergence results in the limit of infinitely many particles: The 
discrete action functionals T-converge to a continuum action functional acting 
on probability measures of particle trajectories and discrete action minimizers 
converge to action minimizers of the limiting continuum functional. In view of 
our numerical convergence results in Section [5] for systems as in [8] , we note 
that our assumptions allow for maximal flexibility in sampling approximating 
marginal measures as we do not assume that these discrete measures are derived 
by some particular coarse graining procedure. 

As an application, Section f3 . 2 1 provides an analysis of the limiting continuum 
mass transportation problem with Lagrangian cost: Minimizers of the limiting 
functional and, more generally, limiting distributions of stationary points are 
investigated and shown to be concentrated on minimizing orbits of the Euler- 
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Lagrange flow. Section [33] then extends the convergence properties from action 
minimizers to general stationary points of the action. 

In Section U the T-convergence result is extended to systems with un- 
bounded potential energies for small time intervals. Under suitable conditions, 
convergence of stationary points of the action is proved even on long time in- 
tervals. 

Section [5] finally gives a rigorous convergence analysis of the aforementioned 
optimal transportation meshfree method, where the discrete trajectories are 
computed through a numerical quadrature formula (the basic midpoint rule). 
This is achieved by extending the previous results to a time discretized set-up 
and proving convergence when simultaneously the time step converges to and 
the number of particles diverges to oo. 

2 Lagrangian systems 

This section serves to collect some notation and well-known basics on La- 
grangian systems. Proofs of these results can be found, e.g., in [7] or [12]. 

2.1 Single particle systems 

Let (M,g) be a connected complete Riemannian manifold of class C°°. By 
™ = Ugm T ^M, T*M = \J xeM T*M, ttm : TM -> M and n* M : T*M -> M 
we denote the tangent and cotangent bundle and their canonical projections 
onto M, respectively. We assume that there is a (time independent) Lagrangian 
L : TM — > R, (x,v) t-t L(x,v), satisfying the following set of classical condi- 
tions: 

(i) Smoothness: L G C°°{TM ; R), 

(ii) Strict convexity: For each compact set K C M there exists some positive 
constant c$ such that V^,L(x, •) > cog x for all x G K. 

(iii) Boundedness: There are constants c±,C2 > such that L(x,v) > c\\v\ — C2 
for all (x, v) G TM. 

By © , L is uniformly strictly convex in v uniformly on compact subsets of M 
and so hin^^oo L ^ v '^ = oo for every x G M uniformly on compact subsets of 
M. Thus, in particular, L is a Tonelli-Lagrangian (see, e.g., [7]). 

In fact, one could more generally also consider explicitly time dependent 
Lagrangians L : [a, b] x TM — > R. Under a suitable additional assumption (see 
below) the necessary modifications are straightforward so that we will not dwell 
on this point. On the other hand, Assumption (iii) is too restrictive for certain 
applications. For this reason we will in Section U] also investigate systems which 
instead of (ii) and (iii) satisfy 

(ii') Uniform strict convexity: There exists some positive constant cq such that 
V^L(x, •) > c g x for all x £ M. 
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(iii') Growth condition: There are constants ci,C2 > and a reference point 
xq E M such that L(x,v) > ci|t>| 2 — 02(1 + d 2 M {x, xo)) for all (x,v) £ TM, 
where (1m denotes the (geodesic) distance on M. 

Throughout Sections [2] and [3l however, Assumptions (i), (ii) and (iii) are as- 
sumed to hold. 

Fix an interval [a,b] Cl To 7 6 C ac ([a, b]; M), i.e., absolutely continuous 
curves 7 : [a, b] — > M, we associate the action 

.4(7) := f L^(t),j(t))dt. 

J a 

By Assumption (iii) this integral exists in (—00, +00]. As it will be convenient 
in the sequel, we extend 4 to C([a,b]; M) simply by setting 4(7) = +00 if 
7 : [a, b] — >■ M is continuous but not absolutely continuous. 

Critical points of this action functional with respect to variations that keep 
the endpoints fixed are called extremals. In particular, for every x,y S M there 
exists a minimizing extremal connecting x and y, i.e., a minimizer of 4 among 
all absolutely continuous curves 7 : [a, b] — > M with 7(a) = x and 7(6) = y. 
The value of the action of such a minimizing extremal will be denoted c(x,y). 
Then c : M x M — > R is continuous. By our hypotheses, minimizing extremals 
7 are C°° and satisfy the Euler-Lagrange equation 

In fact, there exists a smooth flow 4>\ on TM such that every minimizing 
extremal 7 satisfies (7(i),7(i)) = 4>f_ s ( , y(s),'y(s)) for s,t£ [a, 6], i.e., 7 lies on 
an orbit of §\ . 

Remark. For Lagrangians explicitly depending on time, the existence of such 
a flow, which action minimizers follow, has to be assumed. 

Although we will state our results below in terms of the Euler-Lagrange flow 
(fit, we recall here that the dynamical evolution can alternatively be described 
through the associated Hamiltonian and its flow. The Legendre transform in v 
of the Lagrangian L defines the Hamiltonian H £ C°°(T*M]R), 

H(x,p) := sup (p ■ v — L(x, v )) . 

veT x M 

By strict convexity, the right hand side is maximized at p = V v L(x, v) and 
one defines the global Legendre transform as the diffeomorphism V V L : TM — > 
T*M, (x,v) 1 y (x,V v L(x,v)). Via the canonical symplectic structure on T*M, 
one associates a vector field Xjj, which is uniquely determined by 

X H (x,p) = (V p H(x,p),-V x H(x,p)) 

in local coordinates for M. The smooth flow associated to Xh is denoted 4>f. 
The global Legendre transform conjugates 4> L and : 

<f)t = {VvL)- 1 o 4>f o V V L. 
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Example. The most prominent example which we will also discuss from a 
numerical point of view in Section [5] is given by 

L(x,v) = -^-\ v \ 2 ~ V(x) 

for some constant m > and potential V E C°°(IR n ;M), which is bounded on 
M = W 1 . (See Section [J] for unbounded V.) The corresponding Euler-Lagrange 
equation is Newton's equation of motion 

mj(t) = -V%(i)). 



2.2 Many particle systems 

We proceed to give an elementary account on Lagrangian systems with finitely 
many particles. Suppose L is a Lagrangian on the manifold M as above. We 
consider a system of N identical particles moving along curves 1 1— > 7»(t) on M 
with initial positions 7i(a) = (with X{ ^ Xj for i ^ j) at time t = a and final 
positions ji(b) = i = 1, . . . ,N. The associated Lagrangian action is 



Av( 7 ) : = 

i=l 

acting on continuous curves 7 = (71, . . . , 7 jv) G C([a, M w ) and so, in partic- 
ular, 

iV 



8=1 ^ a 



for absolutely continuous curves 7 = (71, . . . , 7 jv) G C ac ([a, b]; M N ). The factor 
jr is introduced so as to measure the action per particle. It gives the right 
scaling in our convergence analysis to be discussed below. For the analysis of 
the discrete action functional with fixed N it is of course irrelevant. 

Being interested in the behavior of the particle density, we define the prob- 
ability measures 



N N 

i=l i=l 

and - assuming that the particles are indistinguishable - are led to minimizing 
the action functional An(j) among all 7 6 C([a,b];M N ) such that 

1 N 1 N 

k E 5 n w = ^ N) and ]v S s n w = ^ • 

i=l i=l 
This in turn amounts to minimizing i over the set 

{7 6 C([a,&];M") : 7i (o) 7i (6) = T{xi) Vi 

for some T with T#(i^ = fJ,^ N ' \, 
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where T#/j, denotes the push forward of a measure \i under the mapping T. 

By first fixing T, we see that the optimal 73 need to be minimizing extremals 
connecting Xj and T(xj). Moreover, the minimal value of the action is 

min | jf E c (^> T (^)) : T *^ N) = »l N) } 

= mmy^c(x,T(x))»i N \dx) : T#»W = , 

where c(x, y) is the minimal action of curves connecting x and y as defined 
above. This shows that, corresponding to our initial and final distributions, 
T is a solution the Monge-Kantorovich mass transportation problem with La- 
grangian cost c. Of course, this is a finite problem for T and so an optimal T 
always exists, showing that with 73 a minimizing extremal connecting Xj and 
T(xi), 7 is in fact a minimizer for An- 

If 7 = (71, . . . , 7^) is a minimizer and so 7, a minimizing extremal for every 
i, we set Vi := 7«(a). This way we can write 

(7<(*)>7i(<)) = $-a{ x ii v i) 

and thus 

7i = ir M (j) L - a (xi,Vi). 

Here ^M4> L - a '■ TM — >■ C([a, 6];M) is the (continuous) Euler-Lagrange flow 
mapping 

(X, U) !->• (t (->• TT M 4>t-a( X i V ))- 

3 Convergence to the infinite particle system 

In this section we first prove our main convergence result on the passage to a 
continuum system in the limit of infinitely many particles. We then analyze the 
limiting problem by deriving analogous properties to the finite system with the 
help of our convergence scheme. Finally we investigate convergence properties 
of stationary points. 

3.1 Variational convergence to a continuum system 

In order to pass to a continuum system, we consider the general action func- 
tional 

A(vr) = J A(i) dvr(7) 

acting on the space of Borel probability measures V(C([a, b}; M)). Here the 
space of curves C([a, b];M) is endowed with the metric of uniform convergence 

dr(l,l')= sup d M (l(t),j\t)), 

te[a,b] 

where du denotes the (geodesic) distance on M. (C([a, b]; M), dr) then is a 
Polish space. 
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Note that with the help of Lemmas 13.41 and 13.51 below the direct method 
of the calculus of variations can be applied to show existence of minimizers on 
weakly closed sets of probability measures. In this section we aim at proving a 
general T-convergence result which in addition shows that any minimizer can be 
approximated by almost minimizers of finite systems with pre-assigned initial 
and final distributions. This not only is the starting point for our numerical 
convergence analysis in Section [5j but will also allow for an investigation of the 
minimizers of A by recourse to an analysis of the minimizers of An- 

For any initial and final distribution fj, a and respectively, we define the 
restricted functionals A AlaiAlfc by 




if Pr a #7r = ^ and pr fe #vr = fj, b , 
otherwise. 



Here pr t : C([o, b]; M) — > M denotes the evaluation at time t: pr 4 (7) = 7(t). 
Remark. If [i a N ^ and /ul are as in JI]), then any tt with A (n) (n)(tt) < oo 

Ma ;Mf, 

satisfies pr a 7r = /j,^ = -h YliLi $xi with Xj 7^ Xj for i 7^ j. The proof of 
Theorem 13. ltfiT|) will show that all the results in this paper remain true if we 
impose the additional constraint that A (jv) (n) (vr) = 00 unless tt is of the form 

Ma 1Mb 

V 



i=l 



for suitable curves 7^ with ji(a) = Xi, so that any such tt satisfies 

A „w „w( 7r ) =-4jv(71!---j7jv)- 

Ma )M(, 

Statement of the main results 

Fix two compactly supported absolutely continuous probability measures \i a 
and yn, on M. Let -Ua be as in ([T]) such that there exists a compact set 

in M supporting all these discrete measures. 

Theorem 3.1. If /xi^ — fj, a and /xP^ — v [i^, then A (jv) (jv) Y -converges to 

Ma )Mb 

^Ma.Mb on 7'(C([a, 6]; Af)) wniTi respect to weak convergence, i.e.: 
(i) liminf -inequality: Whenever tt^ tt in V(C ([a, 6]; M)), rj/ien 



liminf A ( iv) w (vr (Ar) ) > A^ ai 

N— >oo Ma iM(, 



and 



(ii) recovery sequences: for any tt £ "P(C([a, 6]; M)) there exists a sequence 
7r^ -^w in V(C([a,b];M)) such that 



lim A ( jv) (^(vrW) = A MaiMb (7r) 

AT— >0O Ma iMf, 



Moreover, can be chosen as tt^ n ^ = ^ W / or suitable curves 

7 f° eC oc ([a,6];M). 



7 



In view of our numerical convergence results in Section [5] for optimal trans- 
portation meshfree methods as developed in [8] , we note that our assumptions in 
Theorem 13.11 allow for maximal flexibility in sampling approximating marginal 
measures: With the exception of requiring a common compact support, the 
only assumption is that /i a - 1 Mo an d /4 — M6- m particular, it is not 
necessary to assume that these discrete measures are derived from the limiting 
continuum measures by some particular coarse graining procedure. 

This T-convergence theorem is complemented by the following compactness 
result. 

Theorem 3.2. Ifir( N > is a sequence of probability measures such that the mea- 
sures pT a #Tr( N ^ are supported on a common compact set and the corresponding 
sequence A(tt^ n ^) of actions is bounded, then ttW is relatively compact with 
respect to weak convergence. 

A standard argument in the theory of T-convergence now implies the con- 
vergence of minimizers. 



Corollary 3.3. Under the assumptions of Theorem \3.1\ 

Jim minA (iv, (n)(tt) = min p (n). 

Moreover, ifir( N > is a sequence of almost minimizers of A ( N ) (n), i.e., 

Ma 

lim A (iv, ( iv)(7r (Ar) ) = min A Ma M (vr) , 

TV— >00 Ma >Mf, 77 

then admits a weakly convergent subsequence. In fact, any weakly conver- 
gent subsequence converges weakly to a minimizer of Aa ai u h - 

Proofs of Theorems [371] and [3721 

In the proof of Theorem l3.1l it will be advantageous to use a metric that metrizes 
weak convergence on V(C([a, b]; M)). We therefore recall that weak convergence 
of probability measures (as in any Polish space) is equivalent to convergence in 
the bounded Lipschitz distance 

^bl(^, v) := sup jy ipd(n-v) : \\<p\\u P < l| , 

where for any Lipschitz continuous ip : C([a,b]; M) — > M, 

IH| u .==^|y(T)l + -P M T ) ry l . 

7 77 ^ d r (7,7) 

We start by proving two (well known) preparatory lemmas on coercivity 
and lower semicontinuity. 

Lemma 3.4. For every compact set K C M and any C > the set 

{ 7 eC([a,b];M): 7 (a)eK, A{i) < C} 
is relatively compact in C([a,b]; M). 
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Note that this notion of coerciveness is stronger than the mild coercivity as- 
sumption asking that the set of minimizing extremals between two compact sets 
in M be compact in C([a, b]; M). As this Lemma is crucial in our compactness 
results, we include the straightforward proof. 

Proof. If 7 € C([a,b\; M) satisfies A("y) < C, then it is absolutely continuous 
with 

C> f L(-y(t),i(t))dt> Cl f \i(t)\dt-c 2 (b-a). 

J a J a 

But then (with the convention of denoting generic constants of different value 
by the same symbol C) 

d M (l(t)^(a)) < C. 

So by the Hopf-Rinow theorem there is a compact set K' C M such that 
7 (i) G K' for all t G [a, b], whenever 7(a) G K and A(j) < C. 

By our assumption on L there are constants c, d > such that 

L(x,u) > c\v\ 2 - d VxeK',veT x M. 

So A("f) < C implies that for a < t\ < t 2 < b 

rti / rti \ 2 

|2 



< [ 2 \j(t)\dt< ( j" 2 \i(t)\ 
Jti \Jti 



dM(7(*i),7(*2)) < / l7(*)|dt< / |7(0I * (t2-h)* <CVh=Ti 



We have thus shown that {7 G C([a,b];M) : 7(a) G K, A{j) < C} is 
pointwise compact and equicontinuous. The Arzela-Ascoli theorem then yields 
that this set is relatively compact. □ 

Lemma 3.5. A is lower semicontinuous on (C([a, &]; M), dr). 

Proof. As 7^ — > 7 implies that Tfc I [ti ,t 2 ] ~~ ^ 7l[ti,*2] uniformly on every subinterval 
[£11*2] with a < t\ < t 2 < 6, by splitting the action integral we may without 
loss of generality assume that all 7% and 7 are covered by a single coordinate 
chart. We may also assume that A(jk) 1S bounded. Then in local coordinates, 

C> f L( 7 (t),j(t))dt>c f b \j(t)\ 2 dt-d 

J a J a 

for some constants c, d > 0, so that 7^ is bounded uniformly in W 1,2 {(a, b); R n ), 
where n = dimM. But then 7^ — 7 in W 1 ' 2 ((a,b);W n ) and by convexity of L 
in v we obtain 

liminf / L( lk (t),%(t))dt> L( 7 (t), 7 (i))dt. □ 

To abbreviate we write A for A„ o „. and A at for A (jv) (jv) in the following 
proof. 
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Proof of Theorem \3.1\ (0 In order to prove the lim inf-inequality, we assume 
that and, without loss of generality, Aj\[(tt^ n ^) < oo, so that, in 

particular, pr t #7r^ = /4 for i G {a, 6} and all A" G N. Thus, also pr t #7r = /Ltj 
for i G {a, 6}. As .4 is lower semicontinuous by Lemma 13.51 the claim now 
follows from the portmanteau theorem: 

liminf f Ah)dTT^h) > [ AMchrM. 

N^oo J J 

(jirj) It remains to construct a recovery sequence for given tt G V(C([a, b]; M)). 
Note that in view of (0), it suffices to show that lim sup^^^ An(tt^) < A(tt). 
In fact, a diagonal sequence argument shows that it is enough to prove that, 
for any e > there exists a sequence ir^ N > with 

limsupd B L(vr (7V) ,7r) < 0{e) and limsup A N {^) < A(ir) + 0(e). 

TV— >oo N^oo 

Let e > and choose a compact set AT C C([a, b];M) with 7r(AT) > 1 — e. 
Then choose disjoint measurable sets V\, . . . , V m of diameter at most e covering 
K such that fit(dpr t Vi) = for i = 1, . . . , m and t G {a, 6}. (Such sets can be 
constructed by first covering K with m balls B £ u{^/i) of radius |, then choosing 
£ i ^ (f ! |) sucn that Va(dB £i (yi(a))) = n b (dB £i (ji(b))) = and setting V\ = 
B £l (71), Vj = -B £l (7t) \ (Vj. U . . . U Vi-i) for i > 2.) If necessary splitting the sets 
Vi, we may furthermore assume that, for t G {a, b}, 

pr t (^) n pr t (^) = or pr^) = pr t (V,) 

for all i,j. Now choose 7i G ^ with A(ji) < inf 7 gVj -4(7) + £■ 
As /it(dpic t Vi) = for t G {a, 6}, we have 

lim ^fV^)) = Mpr^)) 

N— too 

for i G {a,b}, i.e., the number of points in {x\, . . . ,xn} that lie in pr a (V^) as 
well as the number of points in {y\, . . . ,un} (counted with multiplicities) that 
lie in pr b (Vi) is 

iV M W(pr t (^))>iV Y, ^(Vj) + o(N) 

for t = a or t = b, respectively. So it is possible to find disjoint sets W< C 
{1, . . . ,xn} and a mapping T : U2=i ~~ > {yi; • • • >y/v}> which is injective if 
{y%, . . . , y/v} is viewed as a multiset, such that 

WiCprjy*), T(Wi)Cw b {Vi) and = ATtt^) + o(N). 

Now for each x £ Wi choose a curve 7^ with A(^ x ) < A{^fi) + Ce by setting 



Ixii) = 7j 



(&-a)i- (a + 6)e 
6 - a - 2e 
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for t £ [a + e, b — e] and connecting x to 7j(a) on [a, a + e] and 7«(6) to T(x) on 
[6 — e, 6] on geodesies with velocities bounded uniformly in e. For x ^ U2=i W» 
choose T(x) G {yi, . . . , ?/7v} such that T becomes a bijection when {yi, . . . , y^} 
is viewed as a multiset and let j x be any curve connecting x to T(x) such that 
-4.(7a;) < C f° r some constant C. (This is possible since by assumption all initial 
and end points are contained in a common compact set.) 
We define ttW by 



7T 



(AO = ±_ 



1 " 



Clearly, pr t #7r^ = for i G {a,fr}. In order to calculate dBL(7r^j 7r), let 
9? : C([a,6];M) ->• R Lipschitz with ||<p||Li P < 1. Then, with V := C([a,b];M)\ 
\JZi Vt and W = {x 1 ,...,x N }\ \JT=i W>, 

/ ^tt^ - tt) = £ [ 1 X>( 7 x) - f <pdir). 
Now for i = and sufficiently large TV we have 



< ^ + 7r(Vb) < MVo) + e<3e. 



For i > 1 we obtain from dp (7, 7) < Ce for any 7, 7 G {7^ : x G Wj} U 



< 



N 



V(7i) - t(^)v(7») 



AT 



< Ce 



for large N. Summing over all i we have shown that indeed 

lim sup c?bl (tt^^ , 7r) < Ce. 

Finally we have to estimate the value of the action of . 

Aiv(TrW) = / Afr)d*W = E^E ^ 

i=o 

1 m 1 



< c 



#W 



TV 



E 



TV ^ N \jeVi 

i=l 



inf ^(7) + Ce 



< o( i) + 0(£ ) + £(1 + (1))tt(^) inf ^(7) 



i=i 



< y ^(7) dTr + o(l) + 0(e) = A(vr) + o(l) + O(e). 
This concludes the proof. 



□ 
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We proceed to prove Theorem [37 

Proof of Theorem\3M Suppose A(tt^) = J _A( 7 ) ck^ < C. By Prohorov's 
theorem it suffices to show that the sequence ir^ N ' is tight. Let e > 0. As by 
Assumption (hi) A is bounded from below by —02(6 — a), we have 



TT 



for otherwise 



W ({ 7 E C([a,6];M) : ^( 7 ) > ° + C f ~ ^ }) < e , 



.4( 7 ) *rW > C + C2 g (6 = Q) • e - c 2 (b -a) = C. 

Since by assumption there is a compact set K C M such that 7r^-a.e. curve 7 
satisfies 7 (a) E AT, we deduce from Lemma 13.41 that there exists a compact set 
K E E C([a,b];M) with 

7rW(ir e ) > 1-e 

for every N. This concludes the proof. □ 

3.2 Analysis of the limiting continuum system 

In the language of the theory of optimal transport, minimizers of A- u b are 
dynamical optimal transference plans from fj, a to fib with respect to the La- 
grangian cost c. These objects have been investigated intensively over the last 
years, see, e.g., [12] for a recent account and, in particular, [4] for a result on 
compact manifolds analogous to Theorem 13.61 below. 

For K\ , K 2 C M we denote by r^ n ^ 2 the set of minimizing extremals 
starting in K\ and ending in K 2 . Our main result on the behavior of minimizers 
of A is the following. 

Theorem 3.6. If it is a minimizer for A« a n b , then 

(i) tt is supported on F^aK f or = supp(/x a ), Kj, = supp(/i?,) and 

(ii) there exists a measure r\ on TM such that tt = TTM^-aT^V- 
We begin with the following (well known) preparation: 

Lemma 3.7. If K\,K 2 C M are compact, then T™ n K2 is compact in C([a, b]; M). 

Proof. As the continuous function c is bounded on K\ x K 2 , by Lemma 13.41 
it suffices to show that is closed. Suppose T™"^ B 7^ — > 7. Then 

7(a) E K\ and 7(6) E K 2 . If 7 were not a minimizing extremal, there would be 
7 with 7(a) = 7(a), 7(6) = 7(6) and, by Lemma 13.51 



liminf.A(7 fe ) > A{i) > A(j). 

But then there exists an e > such that for all k sufficiently large one can 
construct a curve 7 fc ]£ by setting 

. / (b-a)t-(a + b)e 
7fc ^ )=7 l b-a-2e 
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for t G [a + e, b — e] and connecting 7fe(o) to 7(a) on [a, a + e] and 7(6) to 7fc(6) 
on [6 — e, b] suitably, so that A(^fk,e) < g"^ ) contradicting the minimality 
of 7^. □ 

Lemma 3.8. Let Ki,K 2 C M be compact. Then there exists lo : [0,oo) — > R 
continuous with uj(0) = and uj(s) > for s > swc/i i/ioi 

^(7)>c(x,y)+a;(dist(7,r™V 2 )) 

/or every 7 G C([a, 6]; M) wii/i 7(a) = x G K\ and 7(6) = y G i^- 

Proof. It suffices to show that for every e > there exists 5 > such that 

dist( 7 , r™V 2 ) > e => ^(7) - c(x, y) > 5 

for every 7 G C([a, 6]; M) with 7(a) = x £ K\ and 7(6) = y G i^2- 

If this were not the case, then there would be an e > and sequence 
7 ( fc ) with 7 «(a) = x fe G K u 7 «(o) = y fc G K 2 , distfrW, Fg™^) > e and 
-4(7^) — c(xfc, yfc) — > 0. c is continuous and, in particular, bounded on K\ x K 2 . 
So .4.(7^) is bounded and, by Lemma I5T41 thus has a convergent subsequence 
(not relabeled) 7W — >■ 7 with 7(a) = x := lim^oo x^ G K\ and 7(6) = y : = 
lim^oo^fc G K 2 - 

Now the lower semicontinuity of A implies that 

.4(7) < liminf A{^ k) ) = lim c(x fc ,y fc ) = c(x,y), 

k— >OD 

which shows that 7 is a minimizing extremal connecting x and y. This contra- 
dicts the fact that 

dist( 7 , IX> 2 ) = lim dist( 7 (fc) , rg> a ) > e. □ 

k— »oo 

Proof of Theorem \3. 9{ i) . We again write A for An a u b and A^r for A (jv) (jv). 

Let 7T be a minimizer of A. For each N G N choose point sets {x\, . . . , xn} C K a 
and {yi, . . . , j/at} C iff, such that 

1 N 1 N 

= n ^ Sxi ^ Ma and ^ ] = N^2 6 y^ Mb - 

i=l i=l 

Suppose = ^ Si^i with 7^(0) = Xj and 7i(6) = T(xj) is a recovery 
sequence for it with respect to the T-convergence of A^r to A. With ui as in 
Lemma 13.81 for K\ = K a and K 2 = we obtain (cf. Section \2.2\ for the second 
inequality) 

1 - 

A(vr) = lim A N (nW) = lim j=^A(ji) 

i=l 

1 ^ 

> - lim sup 2 (c(xi,r(xi)) + w(dist( 7i ,r^ JC J)) 

> lim minA^) + lim / W (dist( 7 , T£^J) okW( 7 ) 
= mmA(n) + J W (dist( 7 , rg^J) dvr( 7 ). 
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Thus, dist(7, r^ 1 ^ ) = 7r-a.e. and the claim follows. □ 

For the proof of Theorem 13 .6lf|TT|> we will need some finer estimates on recov- 
ery sequences. 

Lemma 3.9. For any two compact sets K\,Ki C M there exists a constant 
C > such that \ j(t)\ < C for all 7 G and t G [a, 6], 

We include the short proof of this well known estimate. 

Proof. For every 7 G r^^V it follows from 

rb 

L(7(t), 7 (t))dt < max c(x,y)<C 

(x,y)£KixK 2 

and the lower boundedness of L that there is some t 7 with |7(i 7 )| < C. But 
then (7(4), j(t)) = 4>\-t (7(^7)7 7(^7)) i s bounded uniformly in i and 7 since 
(7(t),7(t)) solves the Euler-Lagrange equation. □ 

Lemma 3.10. Suppose n is a miniraizer of A fla!flb and let K a = supp(/i a ) ; 

= supp(/ib)- XTien /or every e > f/iere exists a probability measure n £ = 
Jf Eili &w G 6]; M)) supported 

d B L(vr £ ,7r) < e. 

Proof. Choose a recovery sequence ir( N ^ precisely as in the proof of Theorem 
I3.6lf lj). so that 

Im^j ^(dist(7,r£> b ))dvrW(7) =0. 

Setting G = {7 : dist(7, T^ n ^ ) < e}, for iV sufficiently large we thus have 
ft\ N )(G) > 1 — e. In particular, for every 7$ G supp^^)) D G there exists 
7i G with d r (7i,7i) < £■ For 7, G supp^W) \ G we choose 7, G r™ 1 ^ 

arbitrarily. Accordingly we define = i ^ i=1 % • 

Now for any Lipschitz (/? : C([a, b}; M) — )• M. with ||y||Lip < 1> 

JipdfrW-nW)^}. y, \<pM-vm+^ E (1^(7)1 + 1^(7)1) 

2#{i : 7, £ G} 
< e + 1 ^ - — 1 < 3e. 

This shows that d,Bh(^ N \ tt^) < 3e. But is a recovery sequence for 7r, 
whence also dBh(^ N \^) < e for large iV. It follows that c^bl^^i tt) < 4e. □ 

Proof of Theorem \3. &{ ii). Let 7r be a minimizer of A^ aiW) . By Lemma f3.10l there 
are measures = X^=i 5 7l with 7, G T^ 11 ^ and such that tt^ n ^ — ^ 7r. 
With = 7i(a) and V{ := 7i(a) we can write 

7; = 4> L - a ( x U v i)- 
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So the measures 

1 N 

i=l 

on TM satisfy = itM ( l )L -aH ir l <yN ^ ■ Note that \v%\ < C for some constant 
C independent of N by Lemma 13.91 Being supported on a common compact 
subset of TM, there exists a subsequence (not relabeled) such that rf N ' — 77. 
But then also itM't'—aH 17 !^ ~^ 7r M0^L a # r ? an d hence ir = TTM4> L - a #r]- C 

3.3 Convergence of stationary points 

In this section we investigate the limiting behavior of critical points which do not 
necessarily need to be minimizers. In a sense to be made precise, we will show 
that in the discrete-to-continuum limit extremals converge to a distribution in 
phase space which follows the Euler-Lagrange flow. We consider sequences of 
curves 7^ = (7^ , ...,7^) G C([a, b]; M N ) and their empirical measures 
71"^) = 7^X^=1°" W with uniformly bounded values of the action An (7^) 

which are stationary with respect to variations in 7W. An being the sum of 
single particle contributions, for every i the curve 7$ then satisfies the Euler- 
Lagrange equation 

±(V v L)(ir\t),4 N \t)) = (V,L)( 7 f \t),4 N \t)). 

The following theorem gives a precise version in which sense extremal curves 
converge to a measure supported on extremals that evolves according to the 
Euler-Lagrange flow. For two compact sets K\,K2 C M let Tk 1 ,k 2 De t ne the 
set of extremals 7 starting in K\ and ending in K2 and set K = {7 G 

t Ki ,k 2 ■■ Aft) < C}. 

Theorem 3.11. Suppose 7^ G C([a,b];M N ) is a sequence of stationary 
points for An with An (7^) < C /or some constant C and ^ N \a) G , 
^ N \b) G iff for compact sets K a ,K b C M. Let tt^ = ^ Eili 5 W- ™ e « 

/or a subsequence (not relabeled) for some ir G V(C([a, b\; M)). ir is 

supported on Tx a ,K b and there is a measure 77 on TM such that -k = 'KM(t )L - a H ir l- 

We begin by proving an extension of Lemmas 13.71 and 13.91 

Lemma 3.12. Suppose Ki,K2 C M are compact. There is a compact set 
K C M and a constant c = c(K±, C) such that for all 7 G K2 

7(i) G K and | 7 (t)| < c for all t G [a,b\. 

Furthermore, K2 is closed. 

Proof. By Lemma 13.41 there exists a compact set K C M with j(t) G K for all 
7 G rjjf j^- 2 and t G [a, 6]. The second claim, namely, I7I < c, follows exactly as 
in the proof of Lemma 13.91 
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Now if K 3 7^ — > 7, then on the one hand we have 

Airfi < liminf^l(7 fc ) < C 

by Lemma 13.51 On the other hand, 7 is indeed an extremal: By splitting, if 
necessary, the action functional we may without loss of generality assume that 
all 7^ and 7 are covered by a single coordinate chart. Then in local coordinates 
by the boundedness of |7fc|, in fact also | 7 fc| is bounded uniformly in k since each 
7^ solves the Euler-Lagrange equation and L is strictly convex on compacts by 
Assumption (ii). We therefore have 7^ — 1 7 in W 2,co ((a, b); W 1 ). But then we 
may pass to the limit in the Euler-Lagrange equations for 7^ to obtain that also 
7 solves the Euler-Lagrange equation. □ 

Proof of Theorem \3.11[ By Theorem 13.21 there is a subsequence (not relabeled) 
such that ttW vr for some vr G V{C([a, b];M)). 
Let e > 0. As in the proof of Theorem 13.21 we find 

7r^({ 7 eC([a,b];M):A( 7 )>c e })<e, (2) 

for Ce = c+c 2 (b-a) and go ^(T*^) > 1 - e. Since T% atKh by Lemma GDI 
is closed, we deduce from and the portmanteau theorem that 

tt(T Ki>K2 ) > n(T% ) > limsupvrW^ 2 ) > 1 - e. 

/V->oo 

As e was arbitrary, Tr^r^jCa) = 1- 

Similarly as in the proof of Theorem [3]6]Jii|) with Vi := ji(a) we can write 

for = (71, . . . , 7/y). So the measures 

1 N 

i=l 

on TM satisfy ^1 = vr M Z l a #?? (7V) . By Lemma [3J2] there exists c = c(e) > 
such that \vi | < c for all j with 7$ S T^ ^ 2 , so that 

7? (7V) ({(x,v) eTM :xeK a , \v\ < c}) > ^(T^J > \-e 

by ((21) • Thus being tight, there exists a subsequence (not relabeled) such 
that rf®' — 1 77. But then also TTM^-aifcW ~^ ^M^-a^V an d hence ir = 

KM<p L - a #ri- □ 

4 Systems with unbounded potential energy 

The theory set forth so far accounts for Lagrangians that are bounded from 
below, as required by Assumption (hi) in Section [2j For this reason, we had 
to assume that the potential energy V in the example of page is bounded 
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(at least from above). This is unsatisfactory from an application's point of 
view as not even a simple harmonic oscillator satisfies such an assumption. In 
this section we will show that, under generic assumptions on the Lagrangian L, 
our T-convergence and compactness results extend to systems whose potential 
energy is not necessarily bounded from above, if the time span b — a over which 
trajectories are observed is sufficiently small. Such an assumption is in fact 
necessary as otherwise the action is not bounded from below. The convergence 
of stationary points, however, can still be justified on long time spans under 
reasonable assumptions. 

In the sequel we will assume that L : TM — > R satisfies 

(i) Smoothness: L e C°°(TM;R), 

(ii') Uniform strict convexity: There exists some positive constant cq such that 
V 2 L(x, •) > Cog x for all x £ M. 

(hi') Growth condition: There are constants C\,C2 > such that L(x,v) > 
ci|t>| 2 — C2(l+d| / (x, xo)) for all (x,v) G TM and a reference point xq G M. 

Statement of the main results 

The action A{pf) is defined as before. (Note that by (iii') for 7 G C ac ([a,b]; M) 
the integral j^L{^if),^{t))dt exists in (—00, +00].) Also the many particle 
action An and the general action functional A as well as the initial and final 
measures /j, a N \ fj, a and fJ,h, respectively, are defined as in Section [3j 

Theorem 4.1. Let b — a < \l With the notation of Theorem \3.1[ if (Ma — 1 



4c 2 ' 

fjL a and — 1 fj,b, then A ( N ) (n) V -converges to A„ o „, onV(C([a,b]; M)) with 
respect to weak convergence. 

Again we also have a compactness result: 
Theorem 4.2. Let b — a< \J j^- If^ N ^ is a sequence of probability measures 

such that the measures pr #7r^ are supported on a common compact set and 
the corresponding sequence A(tt^ n ^) of actions is bounded, then ttW is relatively 
compact with respect to weak convergence. 

As a conclusion to Theorems 14.11 and 14.21 we remark that, under the addi- 
tional assumption that b — a < \Jj^, Corollary 13.31 applies verbatim to the 
systems considered in this section. Also Theorems 13.61 and 13.111 remain valid for 
these systems if b — a < \J j^- This follows in a straightforward manner from 
the arguments in Section [S] using that Lemma [4.41 below bounds trajectories by 
their value of the action functional. (Note that the constant is n °t sharp.) 

It is well known that compactness may be lost on long time intervals, ren- 
dering the search for action minimizers meaningless: 

Example. For M = R, L(x,v) = |d 2 — |x 2 with c > m > and a = 0, b = 
y^7r, for every a G R the curve 'jit) = asin(y^t) is an extremal connecting 

2 

to itself with ^.(7) = j(\/mc — c), which diverges to —00 as a — > 00. 



17 



On the other hand, under suitable local bounds on the action, we may still 
formulate a version of Theorem 13.111 for unbounded potentials on long time 
intervals if the Euler-Lagrange flow is complete. 

Theorem 4.3. Suppose that the Euler-Lagrange flow (pt * s complete. Let 
ry(N) £ C([a, 6]; M N ) be a sequence of N-tuples of solutions of the Euler-Lagrange 
equation with ^ N \a) G , 7^(6) G for compact sets K a , Kb C M. As- 
sume that there is a constant C > such that 

A N ^ N \, V] ) = i£ / L^ N \t),^ N \t))dt < C 
i=i Ja 

for some b' > a with b' < min{6, Then = Y^i=i$ W admits a 

V 

subsequence (not relabeled) such that ir( N ^ — 1 tt for some it G V(C([a,b]; M)). 
tt is supported on T%; a x b and there is a measure n on TM such that tt = 

7T M </> L -a#V- 

Proofs of Theorems EE], 0T2] and I3T31 

We begin by estimating intermediate particle positions by the action. 

Lemma 4.4. Let b — a< y^|j ■ F°r every compact set K C M and any C > 

there exists a constant c = c(c\,C2,K,C) such that for all 7 G C ac ([a, b];M) 
with 7(a) G .FT and ^(7) < C 

d M (rf{t),x ) <c Vi G [a, 6]. 

Proof. As 

d 2 M ( 7 (t),Xo)dt < / 24 / ( 7 (a),x ) + 24 / (7(a),7(t))di 



< 2(6 - a)4 7 ( 7 (a),xo) + 2(6 - a) (J | 7 (s)| 

< 2(6 - a)d 2 M { 1 {a), x ) + 2(6 - a) 2 / | 7 (s)| 2 ds, 



by Assumption (hi') on L we have 

rb 



/■D 1-0 

A(j)>a \j(t)\ 2 dt - c 2 (b - a) - c 2 d 2 M (^(t),x )dt 

J a J a 

> ( Cl -2c 2 (6-a) 2 ) / \"f(t)\ 2 dt - c 2 (b - a)(l + 2d 2 M (j(a), x Q )) 

J a 



>y / \j(t)\ 2 dt- V^(l + 4f(7(a),x )). 
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b / r b \\ 



But then for every t £ [a, b] 

d M (7(t) j7 (a))< ^W)\dt<Vb^U \j(t)\ 2 dt 

i 

< (^^^(7) + 2 + d 2 M ( 7 (a), x ))(6 - a)) 2 , 

and so d M (-y(t), x ) < d M (j(a),x ) + d M (l(a),~f(t)) < c for c = c(ci, c 2 , if, C) 
sufficiently large. □ 

Proof of Theorem\4l\ (i) For R > let fl G C°°(M;R + ) with fl (x) = for 
dM(x,xo) < R and 6r(x) > 02(1 + \x\ 2 ) for d^fO^^o) > 21? and define the 
Lagrangian L R (x,v) = L(x,v) + 9r{x). Then L R satisfies Assumptions (i), (ii) 
and (iii) from Section [2j Denote the induced action functionals by A R , A R , 
A Ma,Mb and A %) (Jv)> respectively. 

Ma ;Mfj 

For given k G N, by Lemma 14.41 we may choose R so large that for each 
7 € C([a, b];M) with 7(a) £ UjveN SU PP Mo we have 

^(7) < =>• d^(7(*),x ) < c < i? => ^ R (7) = A(i). 

But then ^4(7) > minjfc, ^.^(7)} for every such 7. 

Now assume that tt^ N ' — 1 7r in 7- > (C([a, 6]; M)) with uniformly bounded 
A^^ (JV) (7r^). Then as in the proof of Theorem 13. 1 1 prj#7r = fi t for t G {a, 6}. 

Ma X 

By Lemma 13.51 ,A is lower semicontinuous and so is 7 H- min{/c, ^4^(7)}, so 
the portmanteau theorem gives 

liminf f A(i) dn {N) (7) > liminf / min{fc, ^(7)} dir {N) (7) 
Af->oo y jV— s-oo J 

> / min{A;,^4 i? (7)}(i7r(7) 



> y min{A;, ^(7)} ^(7). 

As k was arbitrary, by the monotone convergence theorem we obtain 

liminf f ^(7)^^(7) > / Ah)dnh). 
jV^oo J J 

in the limit k — > 00. 

(ii) For the construction of a recovery sequence of a given measure 7r G 
T^CQa, 6]; M)) it suffices to note that all the measures tt^ n ' of the recovery 
sequence obtained in the proof of Theorem 13.11 are supported on curves with 
values in a common compact set, so that with A R as above for R sufficiently 
large 

limsupA { jv) (jv)(vr^) = limsup A^jv) (jv)(tt^) 

jV-s-oo ^ a '^b jV^oo ^ >^b 

< A R a iW (7r) + O(e) = A /Wtb (7r) + O(e) 
with arbitrary e > as in the proof of Theorem 13.11 □ 
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Proof of Theorem \4-S\ This follows along the lines of the proof of Theorem 13.21 
As .4.(7) < C by Lemma H3] implies that d,M(j(t),xo) < c for 7r^-a.e. 7, A by 
Assumption (hi') is bounded from below on the joint support of the ir( N \ As 
in the proof of Theorem 13.21 this implies that for given e > there is a constant 
C > such that 

ttW (L G C([a,b];M) : ^(7) > | Vj < s. 

Now for 7r^-a.e. 7 we have (7(f), a?o) < c(e) if .4,(7) < — again by Lemma 
14.41 So with as in the proof of Theorem 14.11 for R sufficiently large ^.(7) = 
A R ("f) and by Lemma [3741 there is a compact set K £ with 

| 7 E C([M];M) : ^(7) < || C | 7 G C([o,6];M) : ^(7) < || C A £ 

and therefore tt^-PQ > ^^({7 G C([a,6];M) : ^(7) < f }) > 1 - e. The 
sequence is thus tight. □ 



Proof of Theorem \4-^\ Invoking Lemma 14.41 we may argue as in the proof of 



Theorems 13.111 and 14.21 to find some c e > with 

ttW ({7 G C([a,b];M) : A(-y\ [atV] ) > c £ }) < e. (3) 
Similarly as in the proof of Theorem !3.6|[fTT|) setting Vi := ji(a) so that 

for 7^' = (71, . . . ,7iv), the measures 

N 

A 



' AT 



i=l 



on TM satisfy tt^ = -ir M (f) L _ a #V {N) ■ 

Applying Lemmas 13. 121 and 14.41 on the time interval [a, b'] we find c = c(e) > 
such that < c for all i with -4.(7i) < c £ , so that 



r/ w ({(x,v) G TM : x £ K a , \v\ < c}) > 1 - e 

by (J3]). Thus being tight, there exists a subsequence (not relabeled) such 
that rj. But then also TTM^-ajfcV^ ~ 1 ^M^-a^V as measures on 

C([a, 6]; M) and hence 7r = ■KM(t> L - a H ir l- d 



5 Numerical schemes 

In this section we show how our results can be applied to investigate convergence 
properties of optimal transportation meshfree methods as described in [8]. In 
such a scheme also the time is discretized, and we need to extend our analysis 
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to a set-up where simultaneously the time step converges to while the number 
of particles tends to infinity. To this end, by identifying the time discretized 
system as a perturbation of the systems analyzed in Sections [3J and H] we will 
be able to reduce to our earlier results. While such a reduction is possible for 
many numerical quadrature schemes, by way of example we concentrate on the 
midpoint rule here. For the sake of simplicity we also restrict our analysis to 
Lagrangians L : IR n x R n — > R of the form 

L(x,v) = -^-\v\ 2 — V(x) 

with m > 0, where V G C°°(M n ) satisfies |V(a;)| < c%{\ + |x| 2 ) for some C2 > 0. 
Note that by conservation of the energy y|7(*)| 2 + V(j(t)) for solutions 7 of 
the Euler-Lagrange equation, in particular one has |7(£)| < C(l + |7(i)|) so that 
by Gronwall's inequality the associated Euler-Lagrange flow is complete. 

For a given triangulation Th = {(t 3 -_i, Tj) : j = 1, . . . , 1} of the time interval 
(a, 6), where a = tq < t\ < . . . < t\ = b with \tj — r 3 -_j| < h, the corresponding 
action A shall be discretized by the midpoint rule as 

^]( 7 ) = £ - l7(Tj) 7 (Tj ' l)|2 - v ( jM±^A) (t . _ T ._ l)t 

j = i ^ T j T j-i v 2 y 

In fact, A [h] is a variational integrator only depending on the nodal values j(rj). 
As it is convenient, we may and will assume that A^^) is finite only for 7 
piecewise afflne subordinate to Th, in which case we may also write 



^'(7) = [ f \m 2 dt-jy^ (^^)) (r 3 - w 



Statement of the main results 

Suppose n a , /ift, fij^ ,1^1* are as in Section [3] and let A^, Aj£j )W) and A 

denote the time discretized action functionals on V(C([a, b]; R n )) obtained by 
replacing A by A^ . 



(JV) (JV) 



Theorem 5.1. Suppose V is bounded or b — a < y 3^-- Let hjy > with /i/v — >■ 

as N — > 00. With the notation of Theorem \3.1[ if /ui — 1 [A a and fj,b, 
then A^l (N) Y -converges to Aa «, on V(C([a,b];M. n )) with respect to weak 
convergence. 

Again we also have a compactness result: 

Theorem 5.2. Suppose V is bounded or b — a < \ Let hjy > loit/i 

/i 7v — > as N — > 00. If tt( N ' is a sequence of probability measures such that 
the measures pr a #Tr( N ' are supported on a common compact set and the cor- 
responding sequence A^l (ir^ N ' ) of actions is bounded, then ttW is relatively 
compact with respect to weak convergence. 
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Again we note that, as a conclusion to Theorems I5.1l and l5.2l if V is bounded 

IN) (?• 

Ma 



or 6 — a < . I g^-, then Corollary 13,31 remains true if A ^jvj ^jv) is replaced by 



tm • For a result on stationary points corresponding to Theorem I3.1H 

Ma ',(4 

however, we do have to argue differently, as the discretized trajectories only 
solve a discretized version of the Euler-Lagrange equations and, in particular, 
are not extremals of A. (We also remark that the constant is n °t sharp.) 

Theorem 5.3. Let 7W G C([a,b}; (M n ) N ) be a sequence of stationary points 
for A^ N ^ with 7^'(a) G for a compact set K a C W 1 , where > with 
h]\r — > as N — > oo. Assume that 

^ Nl (7 (JV) |[a, 6 '])<C 

/or some constant C and b' > a with b' — a < y^gf^j- Then there exists a 

subsequence of = j?Yli=i$ W (not relabeled) such that it for 

some it G V(C([a,b];M. n )). it is supported on Tx a ,K b and there is a measure r\ 
on M. 2n such that it = TT^^-a^V- 

Proofs of Theorems EE EJ and RT5I 

Similarly as in Lemma 14.41 we have: 



Lemma 5.4. IfV is bounded orifb—a< y then there exists a constant 
c = c(m,C2,R) such that for all 7 G C ac ([a, &]; R n ) toitt (7(a) | < i? 

| 7 (t)| 2 dt <c(l+„4 [/l] ( 7 )) and |7(t)| <^ W (7) + c Vte[a,6]. 



Proof. If V is bounded by |V| < C, say, we clearly have 



In case V is unbounded, we define 7 : [a, 6] — > 1" by 7(0) = 7(a), 

l(Tj)=l(Tj), 7 Kl ' =7 H 5 ' =7' 



for j = 1, . . . , I and affine interpolation. As in Lemma 14.41 we then find 
l7(t)| 2 dt < 2(6 - a)R 2 + 2(6 - a) 2 / | 7 (s)| 2 ds. 

J a 

Now by construction 7 satisfies 

£ v (t ( a± ^)) - 1-1) = ^ E £ST mo) * 

<2c 2 / (l + |7(t)| 2 )dt 

«/ a 
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as well as 



\j(t)\ 2 dt = - / m)\ 2 dt 



So 



"m 

> I - - 4c 2 



/ \m\ 2 dt 

J a 

) J \^(t)\ 2 dt-2c 2 (b-a)(l + 2R 2 



fi(t)\ 2 dt - 2c 2 {b - a) -2c 2 
b 



> 



in 



\^{t)\ 2 dt - ^/^{l + R 2 ) 



^ j\^t)\ 2 dt-^- 2 (l + R 2 ). 



This proves the first estimate. The second estimate then follows, possibly 
after enlarging c, precisely as at the end of the proof of Lemma 14.41 □ 



The following lemma provides a well-known basic error estimate for piece- 
wise affine interpolations. 

Lemma 5.5. For every piecewise affine 7 subordinate to Th, 

l 



V(i{t))dt-Y,V\~ 



T 3 + T j-1 



T 3 ~ T 3-h 



<h 2 sup \V 2 V(-f(t))\ / \j(t)\ 2 dt. 

a<t<b J a 

Proof. This immediately follows from a Taylor expansion on every subinterval 
(Tj_i,Tj), on which 7 is affine: 



V( 1 {t))dt-V[ 1 



r 7-l + T i 



( r i - t j-i) 



< 



\ T j - T i-i\ 
24 



sup | 7 T (t)W( 7 (t))7WI 



<h 2 sup |VY( 7 (i))| \i{t)Ydt. 

JTj-l 



□ 



Proof of Theorem \5.1[ (i) We define A R as in the proof of Theorem 14.11 and for 
given k > choose R so large that by Lemma 15. 



^[M( 7 ) < k | 7 (t)| <k + c<R ^(7) = ^(7). 
By Lemmas 15.41 and 15.51 we obtain that for those 7 in addition 
|^[M( 7 ) _^( 7 )| < ^ 2 C (1 + k) sup |V 2 F(x)|. 

|x|<fe+c 
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But then, if converges to tt in V(C([a,b]; M)) weakly with uniformly 

bounded action, so that pr t #7r = fj^ for t G {a, &}, 

X^Ct)^ > J mm{k,A [hN] {-f)}d7r N > J (mm{k,A R (j)}-Ch%) dv N 

for a constant C = C(k). Letting first N — > oo, — > and then k — > oo, by 
the portmanteau theorem and monotone convergence we thus arrive at 

liminf / A [hN] h)dn N > [ Ah) dir 

N^oo J J 

as in the proof of Theorem 14.11 

(ii) In order to provide a recovery sequence for tt £ V(C([a, b\; M n )) we only 
need to observe that replacing the curves j x constructed in the proof of Theorem 
13.11 with their piecewise afline interpolation subordinate to Th N only introduces 
negligible error terms in d-Bhi^^h tt) and in A(7r^) in the limit — > 0. Still 
being supported on a common compact set, this defines a recovery sequence 
also for unbounded V, as shown in the proof of Theorem 14.11 □ 

Proof of Theorem\EM As A^h) < C by Lemma IO implies that \j(t)\ < 
C + c for vr^-a.e. 7, A^ hN ^ is bounded from below on the joint support of the 
it( n \ As in the proofs of Theorems 13.21 and 14.21 we therefore have 

tt^ ^| 7G C([a,6];M n ) :^]( 7 ) > j^j < e. 

for some constant C > 0. As at the beginning of the proof of Theorem 15.11 
we see that Lemma I5H implies \A^h) - ^(7)) < 1 for A^(j) < ^ and 
sufficiently small h and large R, so that Lemma 13.41 yields a compact set K e 
with 

1 7 e C([a,b];M. n ) : ^(7) < j j C | 7 6 C([a, b}; R n ) : ^(7) < j + l| 

CK £ 

and^ N \K £ ) > fW({ 7 6 C([a, b}; R) : A^h) < f}) > 1 - e. □ 

We finally turn to the proof of Theorem 15.31 If 7 is a stationary point of 
A^- h \ then it satisfies the discrete Euler-Lagrange equations 

m 7(Tj) -7(tj-i) _ m 7( r i+i) - Tfa) 



73 - 1 Tj+i - Tj 

T i- T i- l vv ( 7(^-i) + 7(^0 ^ + ^+1 ~ 7j w f 7(Tj) + 7(Ti+i) 
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Lemma 5.6. Suppose 7 is a stationary point of A^ with \~f(t)\ < R for all 
t G [a, b] and A^ < C . Let 7 be the extremal of A with 7(a) = 7(a) and 
7(a) = ^( a+/ 7 ^ . T/ien i/iere exists a constant c = c(R,C) such that 

h(t) - j(t)\ + \i(t) - *f(t)\ < ch 

for all t £ [a, b]. 

Proof. If t G (rj-x,Tj), we have 

j, t \ _ 7(ti) -7(to) + / 7(Ti+i) ~ Tpg) _ 7_fa) - 7(^-1) 
_7(n)-7(«) ig^/T^t^)^-,) 



ti — a 2m ^— ' V V 2 



+ vr( 7(n) + 7(n- 1 ) )(r _ rj) 



7(Tl) " 7(a) 1 ^\v(j(s))ds + 0(h) 



t\ — a m 



by Lemma 15.51 with V replaced by W. Note that by Lemma 15.41 the error 
term only depends on C and R. Since 7 solves the continuous Euler-Lagrange 
equation, we also have 



1 



7 (*) = 7(a)-- / VV(j(s))ds. 
m 

Then 7(a) = ^( a+ft 7( Q ) an( j ^.( a ) = 7(a) imply 

|7(*) " 7(*)l < - T I W( 7 (*)) - W( 7 (*))| ds + C7i 
m ./„ 



<C / |7(s)-7(s)|ds + C/i 

<C I f \j(r) -^f(r)\drds + Ch 

J a J a 

= C I {t- r)\j(r) - 7(r)| dr + Ch 

J a 

< C f |7(r)-7(r)|dr + Ch 



for some C > and all a < t < t* := min{min{s : |7(s)| > 6}. Gronwall's 

inequality now yields |^y(t) — 7 (i)| < Che c ^~ a ^ and because of 7(a) = 7(a) for 
suitable c > we thus have 

\j{t)-i(t)\ + W)-j(t)\ <ch 

for these t. But then |7(t*)| < | 7 (i*)| + ch < R + l for sufficiently small h and 
so in fact t* = b. □ 
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Proof of Theorem 15.31 As in the proof of Theorem 15.21 we find 

ttW ({7 G C([a,b];R n ) : -A^M) > c £ }) < e (4) 

for given e > and suitable c e > 0. Consider the discrete Euler-Lagrange 
mapping <j).l a , where '• ^ 2n ~~ ^ ^ 2n maps (zc, v) to the solution of the 

discrete Euler-Lagrange equation 7 with 7(a) = x and 7 ^_^ a ^ = u, so that 
with Vi := !M0JzIiW f or j( N ) = (^ _ _ _ ; ^) we can wr ite 

7i = 4> L -a N ( x ii v i) 

and the measures 

1 * 

i=l 

on M 2 ™ satisfy = TTnn(f> L l a N jfcr]^ ■ By Lemmas 13.121 and 15.61 (applied on 
the time interval [a, b']) there exists c = c(e) > such that \vi\ < c for all i with 
«4 [M (7i|[a,6']) > c e , so that 

»7 (JV) ({(x,t;) £l 2n :ie |«| < c}) > 1 - e 

by (HJ). The sequence rj( N > thus being tight, there exists a subsequence (not 
relabeled) such that — ^ r/. Noting that Lemma [5 . 61 implies that <^2„ — >• </>^_ a 
uniformly on compact subsets of ~R 2n , we finally obtain that TTM n <p L l a i^ r l^ N ^ ~^ 
TTR«(f> L -a#'n and hence ir = n^cj) L _ a #r]. □ 
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